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THE INTERACTIONS BETWEEN NITROGEN AND OXYGEN MOLECULES 


By Willard E. Meador, Jr. 


SUMMARY 

A mathematical analysis is given oj the delta 
junction model for atomic interactions for the pur- 
poses of (1) establishing conditions for which the 
procedure is applicable, and (2) obtaining physical 
insight into the reasons why the method , simple 
though it is } yields potential curves which are in such 
good agreement with experimental data . Lippin- 
cott's original model is then extended to include the 
different effects of K- and L -shell electrons in molecu- 
lar bond {or antibond) formation , and a screening- 
dependent parameter is introduced in the expression 
for the delta function strength. In addition , modifi- 
cations are made which make the model more general 
in application — in particular , in regard to hetero- 
nuclear molecules. 

Modified versions of the delta function model , 
together with general valence bond and molecular 
orbital theories and a reasonably extensive treatment 
of resonance , dispersion , and configuration inter- 
action phenomena , are applied to selected excited 
states of the N 2 , NO, and 0 2 molecules . The results , 
in conjunction with known spectroscopic data and! or 
calculations of the ratios of exchange integrals , are 
then used to find curves representing N 2 — N 2 , 
N 2 -0 2 , and 0 2 — 0 2 interactions , chosen because of 
their importance in problems of aerophysics. Al- 
though the absolute accuracy of the latter potentials is 
difficult to ascertain, they are at least consistent with 
available scattering and viscosity measurements. This 
is especially true of the N 2 —N 2 calculation in which 
remarkable agreement with scattering experiments is 
obtained . It is further believed that the points 
covered in this paper will prove useful in future 
investigations of the interactions between ions and 
neutral species. 


INTRODUCTION 

Interactions between oxygen and nitrogen mole- 
cules are of fundamental importance in the study 
of atmospheric transport properties and other 
phenomena associated with man’s venture into 
space (ref. 1). As a first step in the theoretical 
investigation of scattering cross sections, for ex- 
ample, it is necessary to have a fairly detailed 
knowledge of the functional form of the interaction 
potentials between the elements of the gas under 
consideration. However, the standard methods of 
molecular quantum mechanics (e.g., the ITeitler- 
London (ref. 2) and molecular orbital (ref. 3) 
approaches), in spite of the very lengthy mathe- 
matical procedures involved, do not result in 
sufficiently accurate curves for quantitative analy- 
sis. It is therefore necessary to develop new 
techniques for the calculation of those curves or 
else resort to their determination from experi- 
mental scattering and viscosity data. The latter 
can be extremely difficult and, while acceptable 
as far as the end result is concerned, is certainly 
not as satisfying to the theorist as is the former. 

Mason and Vanderslice (ref. 4) have recently 
presented a method for the calculation of inter- 
molecular forces using a one-dimensional model 
in which the nuclear-electronic coulomb potentials 
are replaced by delta functions. Tlius far, most 
of the applications of this model have been made 
on such simple systems as hydrogen (ref. 5) and 
rare gas atoms (ref. 4) because of the spherical 
symmetry and closed electron shells, which imply 
only one possible interaction curve. The results 
are in surprisingly good agreement with experi- 
mental data in view of the simplicity of the as- 
sumptions involved. 
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Despite this success, however, it has been 
generally conceded that the model is really only 
semiempirical and that a cancellation of errors 
plays a large part in the results. The purpose of 
the present paper is fourfold as follows: 

1. To analyze the delta-function model and 
show why it works. 

2. To apply modified versions of the model to 
N 2 —N 2i N 2 -0 2 , and 0 2 — 0 2 interactions. 
The first example is used to iron out many 
of the procedural difficulties and is chosen 
because of the chemical similarity of molec- 
ular nitrogen to the rare gas atoms. Also 
experimental scattering data exists for this 
case. 

3. To investigate the effects of a variable 
screening parameter. 

4. To investigate the importance of disper- 
sion forces, especially in regard to their 
short-range cut-off behavior. 

LIST OF SYMBOLS 

a radius of the outermost electronic shell 

in an atom 

a, A, B, constants which appear in the approxi- 

c,d mate wave functions and interaction 

potentials 

a,b coefficients of linear combination of 

atomic orbitals 

c variable used in the delta function 

model and related to the atomic 
energies at infinite internuclear sep- 
aration 

d molecular bond length 

E energy 

g delta function strength 

H total Hamiltonian operator 

J atomic ionization potential 

I n ionization potential of atomic hydrogen 

J Li ionization potential of atomic lithium 

J tj exchange integrals between atomic 

orbitals i and j on different atoms 
1c ratio of squares of overlap integrals; 

also resonance parameter 
n number of effective electron pairs in a 

diatomic molecule 

N symbol for atomic nitrogen; also an 

atomic orbital centered on nitrogen 
when used with the subscript x, y, 
or z 


O symbol for atomic oxygen; also an 

atomic orbital centered on oxygen 
when used with the subscript x, y , 
or z 

j) exponent defined in connection with 

the delta function strength 
P probability of the resonance state 

0 + — N“ 

r distance of an electron from its own 

nucleus; also interatomic separation 
E distance between centers of mass of 

molecules 

R m position of Van der Waals minimum 

Su overlap integral between atomic or- 

bitals i and j centered on different 
atoms 

V potential of electron; also interaction 

potential for diatomic molecules as 
defined in text 

V 0 constant which appears in the expo- 

nential curve fit of several interaction 
potentials 

x, y,z Cartesian coordinates of an atomic 

electron ; also variables defined in 
connection with the overlap integrals 
Z effective nuclear charge of L-shell 

electrons 

a parameter in the modified Buckingham 

empirical function 
a,j3 spin functions 

parameter in the Hulburt-Hirsclifelder 
empirical function 

7 variable related to the resonance pa- 

rameter 

e half width of square well potential 

(approaches zero to form delta func- 
tion); also depth of Van der Waals 
minimum 

X resonance parameter 

y, v t <p elliptical coordinates 

p distance between delta function centers 

r volume element (dr) 

interaction potential between mole- 

cules 

V interaction potential between mole- 

cules as obtained from an expo- 

nential curve fit 

\p wave function (subscripts o and n 

signify ground and nth excited 

states, respectively) 
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os frequency of cosine wave function (see 

°q- (4b)) 

Energies are given in units of twice the ioni- 
zation potential of atomic hydrogen, or in electron 
volts, as noted. Distances are given in units of 
first Bohr radius of atomic hydrogen, or in ang- 
stroms, as noted. 


THE DELTA-FUNCTION MODEL 


HYDROGEN-LIKE ATOMS 

The wave equation for a ground state li 3 T drogen- 
like atom is, in atomic units, 


where 


-U±( . 

2 r 2 jdr \ 


d 

dr 


;)+2Zr 


=Eo'P„ 





Z>\ 


1/2 


(la) 


(lb) 


and Z is the atomic number of the nucleus. 

Now, considering the application of the Lap- 
lacian operator in Cartesian coordinates to ^ 0 , 
we obtain 


b%, _/Z 2 X 2 Zx 2 
dr 2 ~ \ r 2 r 3 



(2a) 


is the crux of the simplification with the result 
that the wave function of the atom may be 
represented in one dimension b}" the dotted curves 
shown in sketch (a). 



and 


]_ 

2 






It is apparent then that the simple mathematical 
procedure of reducing the problem to one dimen- 
sion by taking the components y and z to be zero 
in yp 0 and equation (2a) results in 



This last equation is obviously equivalent, at 
least insofar as the mathematical eigenvalue 
problem is concerned, to equation (la) but with 
the very important physical exception that a 
potential energy function is not included. The 
fortuitous cancellation of the second and third 
terms on the right-hand side of equation (3a) 


Sketch (a). — The hydrogen atom wave function in the 
delta-function model. 


The correct excited state energies are also 
easily obtainable from equation (3b) by simply 
using ^*=exp(— Z\x\Jn) } where n is the principal 
quantum number; however, the situation tends 
to become more obscure and further removed from 
physical reality because of t lie fact that the 
true wave functions are not simple exponentials. 

Since it seems possible to get along so well 
without worrying about a potential energy 
function in the one-dimonsional Hamiltonian, the 
question might reasonably be asked, at this point, 
as to why Mason and Vandersliec (ref. 4) bother 
about introducing a square well at the position of 
the nucleus and then letting the width of this well 
approach zero and its depth approach infinity 
in such a way as to degenerate into a delta func- 
tion, especially since the correct energy values 
are obtained without doing this and the model is 
but a mathematical construct. The answer is 
that in order to apply the method to the more 
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complicated cases of molecules, it will be neces- 
sary at least that the equations represent some 
sort of physical situation, even though it might 
not be that of an atom, and that they satisfy the 
boundary conditions of general quantum theory. 
It is quite apparent that the latter condition is 
not fulfilled in the above simplified model from 
the mere fact that the slope of the wave 
function is not continuous at the origin, an obser- 
vation which clearly indicates a potential center 
at this point . 

A possible procedure for circumventing the 
above difficulty is to connect the two exponential 
regions (i.e., oxp(cx) for and oxp(— cx) for 

x^>0) with a suitable even function of x, such 
that the slopes are continuous, and then to find 
the corresponding potential from which it can 
originate. Perhaps the simplest such connecting 
function is the cosine, the frequency of which 
may be allowed to approach infinity, and which 
is the wave function of a particle in a one-dimen- 
sional box -lienee the use of the delta function 
type of potential (see sketch (a)). It follows then 
that while the delta functions associated with two 
such atoms cannot, of course, overlap, the corre- 
sponding wave functions can. The situation for 
the individual atom is, in some respects, similar 
to the "tunnel effect” in the a decay of radio- 
active nuclei. 

In other words, consider (he a* axis to be divided 
into three regions with the following set of wave 
functions and potentials: 

exp(c.r), F=0(x< — t) (4a) 

^ n = 7? cos cox, V= — V 0 (—€<x<e) (4b) 

fcn=A oxp( — cr), V=0(x>e) (4c) 

where c is essentially an effective nuclear charge 
and is equal to (— 2Z?) 1/2 in this problem. The 
parameter a? is determined in such a fashion as to 
satisfy the one-dimensional wave equation, in- 
cluding the delta-function potential, and 2e repre- 
sents the width of the square well. Thus 

1 d 2 

— - (cos o>j) + V cos cox=E cos cox (5a) 

or 

u2 =2 (E-V)=2{E+V 0 )=-c 2 +2V 0 (5b) 

The quantity co must be real, of course, otherwise 
would be a hyperbolic function. 


Now V Q may be obtained from the requirement 
that the functions and their first derivatives join 
smoothly at the boundaries of the regions, that is, 


II 


at x= — e 

(6a) 


\f/ U ' — \plll 

at x=e 

(6b) 


This results in w tan coe~ c and, since e may be al- 
lowed to approach zero without loss of generality 
for purposes of simplification, wo have 

a > 2 e = — c 2 e + 21 T 0 t = c (7 a) 

Letting g=c-\-c 2 e, this yields 

F„=f, «»=-c*+2 (7b) 

In this example, < j~c , but this is not necessarily 
true in the case of molecules. This point will be- 
come clearer in a subsequent section of the report 
on H 2 + 

It should be noted here that the shrinking of 
the width (2c) of the potential well to zero auto- 
matically leads to the increase of the depth (y/ 2c) 
to infinity, but in such a way that the product of 
width and depth is equal to g, a finite number. 
Thus, the potential used here is indeed a delta 
function of strength g\ however, it is in no way 
unique, that is, this is not the only way in which 
a one-dimensional analysis can be made physically 
plausible. 

Finally, it must be remembered in what follows 
that, when use is made of this model, no coulombie 
interaction of any kind appears explicitly in the 
Hamiltonian operator of an atomic system. 
Therefore, as follows from the first-order nature 
of this procedure, induction and other second- 
order effects must be added arbitrarily in order 
to obtain a complete description of forces between 
two or more interacting species. Moreover, any 
analysis involving a three-dimensional quantity, 
for example, angular momenta different from zero 
and their interactions with spin resulting in fine 
structure corrections, must be handled in a some- 
what different manner. The loss of this type of 
general physical insight, however, is not an un- 
common occurrence when mathematical simpli- 
fications are made for the purpose of dealing with 
a more specific effect. 
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THE H 2 + AND H 3 MOLECULES 

Following Lippincott’s (ref. 5) method of bring- 
ing two delta-function atoms together to form a 
delta-function molecule, we may divide the x axis 
into five regions with the following wave functions: 



vhi=#i cos ux e<x<— 2 + € ) (8b) 

hn= /^[exp ( — cx) + exp (ex) ] 

(-'+.<*<§-.) («,•> 

^iv=^i 2 cos wx (^~ «<Z<^+«) (8(1) 

tv^Ae-' 1 ( x ~2 +e ) (8e) 

where p is the distance between delta-function 
centers. The positive sign is used in \j/ ul to make 
tlie wave function symmetric in an interchange of 
nuclei, since, by general molecular orbital theory, 
this corresponds to a concentration of probability 
density in the region between nuclei and leads, 
therefore, to the lowest energy. 

Using equation (7b) and applying the appro- 
priate boundary conditions to the region end- 
points results in the equations 

(9a) 

c=(j(\+e~ Cp ) (9b) 

^=-| 2 (l + e - fp ) 2 «-| 2 0 + 2 e -") ( 9 c ) 

Equations (9b) and (9c) may now be solved 
simultaneously to determine c and E as a func- 
tion of the delta-function separation p. The 
parameter (j is, of course, determined simply by 
the requirement that the energy approach the cor- 
rect isolated atom energies as p goes to infinity 
(g= \ for ITU and 2 1/2 for II 2 ). 

Lippincott now makes the further argument 
that the approximation should be considerably 
improved by using “floating” delta functions, 
that is, by allowing the delta function centers to 


be shifted off the nuclear centers. This procedure 
will clearly enable the charge distribution to be 
more concentrated in the region between nuclei 
(p<>= interim clear separation) in the case of 
bond formation and in the regions outside of the 
nuclei (p^>r) in the case of repulsive states. 
This is essentially equivalent to the shifted 
atomic orbital procedure as applied to H 2 by 
Gurnee and Magee (ref. 6) and to TT 3 by Meador 
(ref. 7), and it is also similar to second-order 
perturbation theory in which, for example, hybrid 
Is ami 2 p orbitals are used for IL in the ground 
state (ref. 8). The methods for specifying p as a 
function of r will become apparent in the next 
section. 

THE EXCITED STATE OF N* 

In a recent paper Vandcrslicc, Mason, and 
Lippincott (ref. 9) applied the delta-function 
model to the calculation of the energy of the 
7 2u + state of X 2 , in which all the valence electrons 
are unpaired. The antisymmetric wave functions 
appropriate to repulsive states were used in place 
of the symmetric ones discussed previously and 
resulted in the equations 

c-u( \—e~ Cp ) (10a) 

E=-^{ i-2c-‘0 (10b) 

where n is the number of electron pairs (seven in 
this case). The assumption is made that the 
many electrons may be replaced by a single 
effective charge distribution; more precisely, the 
total interaction energy consists of n times the 
average interaction energy of all pairs of electrons, 
as for t wo hydrogen atoms. 

The interaction energy may now be expressed as 


V(r) = E(r) — E(^) — ng 2 e~ cp 

OD 

where 


P = H-2 ac~ rfa 

(12) 

Tli is last equation is merely a simple 
satisfy the requirements: 

s way to 

(a) p>r (repulsive state), p~^r 

as r— > oo 


(b) p(r=0) = 2a, where a is the average ra- 
dius of the outermost electronic orbit of 
an isolated atom (a— 0.56 A (ref. 10) 
for nitrogen) — the “united atom.” 
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With the idea that the delta-function strength 
(g) may be regarded as proportional to the 
ionization potential of an isolated atom, as implied 
by the earlier scattering investigations of Bloch 
for heavy targets (ref. 11 ), one obtains 

»= 2 < 7 .)« (£)-*»(£) 03 ) 

where / and / H are the ionization potentials of 
atomic nitrogen and hydrogen, respectively. This 
procedure finally yields a potential energy which 
can be represented by the expression 

— 317.8 exp (—2.753 r) (14) 

where r is measured in angstroms, and V is in 
electron volts. 

The present paper is an extension of the above 
work in that allowances are made for the smallness 
of the contributions of the inner or K-shell elec- 
trons to the repulsive forces between atoms (or 
to bond formation in the case of bound states). 
The question is raised also as to whether the 
delta-function strength should be proportional to 
I or to 7 172 , or perhaps to some power m between, 
especially in view of the fact that equation (13) 
signifies a discontinuity in going from hydrogen 
(g^P /2 ) to & more complicated atom As 

a more natural assumption, the introduction of a 
variable parameter in order to make this transition 
less abrupt seems appropriate. In addition, since 
there are only five L-shell electron pairs involved 
in No, it is expected that the final form for g 
should more closely resemble that for hydrogen 
than for the heavy atoms. 

The first extension involves a comparison with 
lithium instead of hydrogen in the equation for 
g . In this case the following system of equations 
results (energy in electron volts, and r in 
angstroms) ; 

p= 0 5292 +2J164<,r/0 ' 56 ( 15a ) 


where p is a number, to be determined, between 
one-half and unity. 

A comparison of this approach with experi- 
mental data will be given below, where it is shown 
that p=% gives very good results. This corre- 
sponds to g= 2 (7/27.206) 1/2 so that the interaction 
energy may be represented (after some manipula- 
tion using Newton's iteration procedure and curve 
fitting) by the expression 

F(r)=253.9 exp(-2.716r) cv (16) 

It should be noted here that the method of approxi- 
mation of perfect pairing, as employed in subse- 
quent sections of this report, involves only 2 p type 
electrons in a description of interatomic forces. 
For this reason, it perhaps would have been more 
consistent to have considered three electron pairs 
in equation (15d), instead of five, and to have used 
boron, instead of lithium, in equation (15b). The 
above treatment emphasizes the shell structure, 
and hence regards only the two 1 $ electrons as 
being Lransnuclear, whereas the latter points out 
the division into subshells and is used only in 
finding relationships between energy states. 

Equations (15) point out clearly the advantages 
of using the delta-function model, at least from the 
standpoint of mathematical simplicity, since the 
wave equation becomes completely separable in 
the electron coordinates, and complicated coulomb 
and exchange integrals, etc., do not enter the 
picture. 

THE N 2 -N 2 INTERACTION 

The approximation of perfect pairing (ref. 12 ) 
describes the interaction between two nitrogen 
molecules as simply the sum of four atomic inter- 
actions. In the case of neutral species (i.e., neg- 
lecting any ionic contributions due to resonance) 
the interaction potential is thus given by a partic- 
ular sum of exchange integrals as follows; 

where 


/ h, y 2 / / V / 5.363 V/*/ 14 
g ' - \27.20(5/ vJ \27.206/ V o.. r 

c=<j(l—e~ Cf ) 

V(r) = 5(27.20 O)^ -<*= 1 36.(tt0*e ~ Cp 


48 V 
363/ 
(15b) 

(15c) 

(15d) 


5 1 the sum over orbitals with paired spins (anti- 

parallel) 

5 2 the sum over orbitals with parallel spins 

5 3 the sum over orbitals with nonpaired or ran- 

dom spins 
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The choice of coefficients, of course, follows 
directly from simple valence bond theory. For 
example, the factor— in the last term is a con- 
sequence of the fact that this represents a com- 
pletely random situation and the degeneracies of 
paired and parallel spin states arc one and three, 
respectively. The exchange energy is, therefore, 

Further assumptions concerning molecular in- 
teractions are: 

1. The exchange integrals arc essentially the 
same as would prevail if the atoms were 
isolated instead of being members of mole- 
cules. The directions of the distortions 
experienced by the electronic charge dis- 
tributions as the molecules approach each 
other certainly tend to make this approxi- 
mation good and, in any case, the error in- 
troduced is probably very small for the 
distances under consideration (2.4 to 
3.2 A). 

2. The resonance contribution is insignificant 
because of the small electron affinity of 
molecular nitrogen. In many respects N 2 
behaves like an inert gas atom so that X 2 ~ 
is extremely unlikely. 

3. Only the valence 2 p x , 2p v , and 2 p z atomic 
orbitals are considered as contributing to 
the exchange forces between atoms. 

4. The coulomb interaction between neutral 
molecules is negligible at fairly large inter- 
nuclear separation in comparison with the 
exchange forces. 

5. Only one N 2 — N 2 interaction curve is con- 
sidered to be of importance, that is, the 
possibility of chemical reaction is ignored. 

The 7 yy state of X 2 may now be represented by 

I (r) ~ Jyy-\~ Jzz) (17) 

and the interaction between nitrogen atoms, either 
or both of which arc members of molecules, by 

V'{r)'*-±(J„+Jn+J») ( 18 ) 

Combining equations (16), (17), and (18), we 
have 

V’(r) F(r)=F,e-*'=127.0e- 2 -™' ev (19) 

The total interaction between two nitrogen 

544424-61 2 


molecules is then given by the sum of the four 
interactions between their constituent atoms, four 
in all, and depends implicitly on their relative 
orientations. For comparison with experimental 
scattering and viscosity data, it is convenient to 
average geometrically (ref. 9) over all possible 
orientations, which yields for the average X 2 — N 2 
interaction 

$ (2? ) — 4 V 0 e ~ hR (¥Rd 2 ) “ 1 [2 ( 6 7?-f 2) (cosh 1 ) 

— 2 bd sinh &d]“^rr ev (20) 

where 

R distance between the centers of mass of 

the molecules 

d bond length of X 2 = 1.094 A (ref. 13) 

— 36.0 //i G London dispersion or second-order per- 
turbation energy (ref. 14) 

Calculated values of from this equation and from 

*'(R) = Ae- nB ~-~ ( 21 ) 

where ^1=658.66 and 72=2.630 for p = 0.5, as 
obtained from curve fitting, are presented in table 
I. Also shown is the energy for £>=0.6 in order 
to show the trend when this parameter is varied. 

Comparisons between the theoretical results of 
this paper and those of Vanderslice, Mason, and 
Lippincott (ref. 9), together with the experimental 
data of Amdur, Mason, and Jordan (scattering 
measurements, ref. 15) and Mason and Rice (vis- 
cosity measurements, ref. 16), are presented in 
table II and sketch (b). The excellent agreement 
would seem to justify the present approach. 

There is, however, a discrepancy at large /? 
where the theoretical curve does not approach 
zero as rapidly as the viscosity data, but this may 
he due partly to the inadequacy of the experi- 
mental results in this region. It is also apparent 


TABLE I. -N a -N a INTERACTION ENERGIES 


7? 

Eq. (20) 
and 
p ^O.G 

Eq. (20) 
and 
p = 0.5 

Eq. (21) 
and 
;j = 0.5 

2. 4 

0. 0272 

1. 0071 

1. 0080 

2. 5 

. 4G37 

. 7720 

. 7723 

2. G 

. 3111 

. 501 1 

. 5005 

2. 7 

. 2404 

. 4508 

. 4507 

2. 8 

. 1812 

. 3429 

. 3432 
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TABLE IT. COMPARISON OF N 2 -N 2 INTER- 
ACTION ENERGIES 



Calculated 

Measured 

R 

Ref. 0 

Eq. (21) 
and 
p =0.5 

Ref. 15 

Ref. 1G 

2. 4 

1. 190 

1. 008 

1. 045 


2. 5 

. 905 

. 772 

. 757 


2. 6 

. G92 

. 591 

. 5G8 


2. 7 

. 527 

. 451 

. 440 


2. 8 

. 400 

. 343 

. 33 G 


2. 0 

. 303 

. 2G1 

. 257 


3. 0 

. 229 

. 198 

. 200 


3. 1 

. 173 

. 149 

. 157 


3. 2 

. 130 

. 112 


0. 093 

3. 3 




. 050 

3. 4 




. 023 

3. 5 




. Oil 

3. 6 




. 005 



o 

/?, A 


Sketch (b) .—Comparisons of N 2 — N 2 interaction energies. 

Unit the two theoretical curves converge at large 
distances because of the decreasing importance of 
distinguishing between Tv- and H-shcll electrons. 
Finally, the fact that the curve of this paper lies 
below the one derived from scattering data for R 
less than about 2.45 is due primarily to the in- 
ability of the model to account for the strong 
repulsive forces encountered when the electronic 
charge distributions effectively overlap; that is, 
the diatomic theoretical curve for the ex- 

cited state does not approach infinity at r equal to 
zero. The specific neglect of coulomb interactions 
in the Hamiltonian operator apparently is valid 
only when the intermolecular separation is large 
enough that the forces (other than exchange) 


between individual particles cancel out. These 
arguments thus make clear the limitations of the 
delta-function simplification and point the way to 
the modifications necessary when one or both of 
the interacting species is an ion. 

Finally, another lower limit 011 the distance 
between molecules, insofar as this theory is con- 
cerned, is the point at which the procedure of 
averaging over orientations becomes invalid. 
Such a point would correspond to a separation 
at least as great as 2 (r//2) = 1.094 A. 

DISPERSION FORCES 

The discrepancy between theory and experiment 
at small R (sketch (b)) may, to a minor extent 
compared to the reasons stated in the preceding 
section, be ascribed to the poor expression for the 
dispersion term. The approximations involved 
in the derivation of this term are such that, in 
short-range interactions, it is usually better to 
neglect dispersion entire]}" in comparison with the 
first-order energies. For this reason it would 
probably have been more realistic to have omitted 
the term — 36.0//? 6 in equations (20) and (21); 
however, because of its not insignificant contri- 
bution to the total energy ( — 36.0 //?* varies from 
— 0.188 at 7i = 2.4 to —0.075 at R = 2.8), an attempt 
must he made fo compensate for the omission of 
the term by allowing the p of equation (15b) to 
be greater than 0.5. 

The best agreement with experiment is obtained , 
and to this extent the entire procedure must be 
regarded as essentially semiempirical, if p is taken 
to be 0.55. The results are shown in table III, 
where column 2 values wore calculated using 
equation (20) without dispersion and column 3 


TABLE ITT, X 2 — X 3 INTERACTION ENERGIES 


* 

Eq. (20) 
without 
disper- 
sion 
and 

p ^ 0.55 

Eq. (22) 

Eqs. (20) 
and 
(15d) 

Eq. (21) 
and 
p = 0.5 

Ref. 15 

2. 4 

1. 000 

1. 003 

0. 989 

1. 008 

1. 045 

2. 5 

. 759 

. 759 

. 755 

. 772 

. 757 

2. G 

. 575 

1 . 575 

. 575 

. 591 

. 568 

2.7 

. 436 

. 435 

. 437 

. 451 

. 440 

2. 8 

. 330 

. 329 

. 332 

. 343 

. 33 G 

2. 9 

. 250 

. 249 

. 251 

. 2G1 

. 257 

3. 0 

. 189 

. 189 

. 190 

. 198 

. 200 
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corresponds to the associated exponential curve 
fit given by 

<£'(/?) = 800.52 exp ( — 2,784/?) (22) 

Column 4 of the same table results from using 
an equation of the form (lod) directly in the 
expression for <£(/?), instead of going through the 
intermediate step of equation (16), and the two 
remaining columns are reprints from table II for 
purposes of comparison. Thus, at least in the 
range 2.5</?<2.8, this last description of the 
interaction does an even better job of matching 
the scattering data, as the average discrepancy 
between theory and experiment at the points 
cited is reduced from 0.014 to 0.004 ev. The 
inadequacy of scattering measurements apparently 
begins to take hold for R greater than approxi- 
mately 2.9. 

In addition, a graphically constructed con- 
necting curve between the present theoretical 
results and those obtained from viscosity data 
could yield a valuable empirical determination of 
the effective reduction factor and short-range 
cutoff of the dispersion forces. This will be shown 
more clearly in connection with X,; — 0 2 . 

For further comparison with the work of Van- 
derslice, Mason, and Lippincott (ref. 9), a cal- 
culation or the interaction energy between nitrogen 
atoms and molecules was made using equations 
(15) in conjunction with j) = 0.55, 

V'(r ) « t V Q^:)=V 0 e~» (cf. eq. (19)) (23) 


TABLE IV.— N-N* INTERACTION ENERGIES 


R 

Eq. (24) 

Ref. 9 

2. 4 

0. 37G 

0. 550 

2. 5 

. 285 

. 418 

2. G 

. 215 

. 318 

2. 7 

. 1G3 

. 242 

2. 8 

. 123 

. 184 

2. 9 

. 092 

. 140 

3. 0 

. 070 

. 107 



Sketch (c). — Comparisons of X — N 2 interaction energies. 


slice, et ah, (ref. 9) curve was calculated from 
$/(7?)=387.8 exp(— 2.733/v) (25) 


and 


$ 1 (/?)=2T ^"“(Wft*)- 1 2(W?+1) sinh (^) 




(24) 


Equation (24), of course, follows from the same 
type of averaging procedure as that used to obtain 
equation (20), and the intermediate stage of equa- 
tion (16) was again omitted; that is, equation (23) 
was used for each value of R but no attempt was 
made to represent the entire range by a single 
exponential function. 

The results are shown in table IV and sketch (c). 
There are, unfortunately, no experimental data 
available to check their validity. The Vander- 


To summarize, the exponential curve fit of 
equation (22) provides a very simple function for 
use as the repulsive part of the potential appearing 
in cross-section and transport integrals. The 
neglect of dispersion forces at large distances, how- 
ever, prevents the usual potential minimum, which 
occurs at 5.160 A as computed from equation (21) 
with the term — 3G.0/7? 6 . The corresponding value 
of the potential at this point is — 0.001 1 cv, whereas 
the Vanderslice, et al. (ref. 9), equations predict 
—0.0001 cv at 5.203 A. On the other hand, in 
high-energy scattering experiments, the effect of 
this attractive contribution is often quite negligible 
and equation (20), without dispersion and without 
further modifications, should give very adequate 
results for this type of measurement. 
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INTERACTIONS BETWEEN NITROGEN 
AND OXYGEN 

METHOD 

For the interaction between nitrogen and oxygen 
atoms, both of which are members of their respec- 
tive X 2 and 0 2 molecules, the approximation of 
perfect pairing yields 

T»=-~ [2 (2) J**+\ Jzz+\ j«+ 2 (2) 

+2^+2 ^ IxJr 2 (•>) 

where the x axis is taken to coincide with the line 
of nuclear centers. Here again only p type 
atomic orbitals are considered and advantage 
has been taken of the small electron affinity of 

N*. 

Since symmetry requirements dictate that J vy 
is identical with J ZZ} the above expression can be 
immediately reduced to 

T T M=-f (J«+2JJ (20) 

In previous papers (refs. 9 and 17) employing 
this approach it has been customary to eliminate 
J xx and Jyy from similar equations by using the 
same approximation of perfect pairing in connec- 
tion with various states of associated diatomic 
molecules, the potential curves of which are 
obtained from spectroscopic data or by adopting 
the delta-function model. Such a procedure was 
particularly simple in the case of N 2 — N 3 because 
of the almost trivial relation expressed in equation 
(19); in other words, only one diatomic state had 
to be considered and this was especially adaptable 
to a delta-function treatment. The present prob- 
lem is considerably more complicated because of 
the following three major factors which enter into 
the nitric oxide calculation: 

1. Any attempt to use a delta-function model 
must take into account the fact that there 
will be two different sets of delta-function 
strengths and “floating” parameters. 

2. The appearance of four p electrons in 
atomic oxygen implies the existence of 
three -electron bonds or antibonds in the 
states of nitric oxide. Thus, a careful 
analysis of the approximation of perfect 
pairing, which works so well in the case of 
N 2 , becomes necessary. 


3. The possible attachment of the “extra” 
oxygen electron to the nitrogen atom to 
form a resonance N~ — 0 + state must also 
he considered. 

The addition of an associated resonance parame- 
ter to the two unknown exchange energies in 
equation (2G) requires a minimum of at least three 
independent relations to solve uniquely the 
0 2 problem. Unfortunately, only the A^IT 
ground state of nitric oxide is sufficiently stable to 
provide enough spectroscopic data for an appli- 
cation of the Rydbcrg-Klein-Rees semi empirical 
method (refs. 9 and 17). 

In a recent paper by Vanderslice, Mason, and 
M aisch (ref. 17) these difficulties arc more or less 
avoided by some arbitrary assumptions concerning 
the resonance phenomenon. In addition, a some- 
what questionable method of obtaining the hound 
4 H state of NO, by using the spectroscopic con- 
stants of the corresponding state of 0+ in a 
Hulburt-TTirschfelder function, is employed. It 
is shown in the present work that a more consistent 
determination of the interaction potentials is 
obtained when a delta-function model is used to 
calculate the resonance parameter. A derived 
relation involving the ratio of exchange integrals, 
together with the aforementioned A" 2 II ground 
state data, is then sufficient to specify completely 
the molecular interaction. 


THE RESONANCE PHENOMENON 

As mentioned above, the ground state A" 2 IT 
potential curve of nitric oxide lias been calculated 
by Vanderslice and co-workers using tin* Rydberg- 
Klein-Rees semi empirical procedure. It lias also 
been shown (ref. 17) that the long-range tail of 
this curve joins smoothly with the following 
Hulhurt-IIirschf elder expression (ref. 18): 

V( 2 n) — G.G09 [(1— c - ^) 2 

+ 0.00780 ^(l+2.G63/3)c- 2 *—l] cv (27) 


where 


0=3.15' 


K- 


— 1 . 1 508 \ 
1.1508 ) 


and r is the distance between atoms (in angstroms). 

If the customary molecular notation (ref. 13) 
is used, the configuration of the seven valence p 
electrons may be described by 

(oOWOOW) 
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where the molecular orbital wave functions arc 
approximately given by 

lb :a0 2 -f 

Ib*:aO a — 6N rj etc. 

and where a and b are coefficients of the linear 
combination of atomic orbitals. The 0 2 and X 2 , 
of course, refer to 2 p z atomic orbitals centered on 
oxygen and nitrogen, respectively. 

An application of the approximation of perfect 
pairing leads to V( 2 U) =J X z-\-J v » phis the contri- 
bution from the remaining lb and IR* molecular 
orbitals. 

This last term, however, concerns three elec- 
trons, and the aforementioned simple theory is not 
adequate to describe such a situation. A very 
feasible extension has been suggested by Linnet t 
(ref. 19) in terms of a mixture of atomic and 
molecular orbitals, which follows directly from the 
properties of matrices and determinants. The 
basis of this procedure is easily seen from a con- 
sideration of the case where one electron is in the 
lb orbital and another is in the lb* orbital and 
where both have their spins in the same direction. 
The wave function for this two-electron problem 
is then given by 

Ib(l)a(l) Ib*(l)a(l) 

lb (2) a (2) Ib*(2)a(2) 

__ [aO a (0+&X a (l)]a(l) [aO a (l)-&N a (l)]«(l) 

\aO t (2) + &N 2 (2) ] a (2) [aO z (2) - iN a (2) ] a (2) 

Adding the second column to the first, dividing 
by 2, subtracting the resulting first column from 
the second, and finally multiplying by 1 fa and 
— 1/ft, gives, apart from a constant factor, 

0 2 (l)tt(l) N a (l)«(l) 

O a (2)a(2) N a (2)a(2) 

In other words, in the description of these two 
electrons, it is irrelevant whether they are said to 
occupy the molecular orbitals or the corresponding 
atomic orbitals, that is, the contribution to V( 2 U) 
is — J vy regardless of the viewpoint taken. 

The remaining electron in the three-electron 
bond must, in order to satisfy the Pauli principle, 
have a spin function j3 and may roughly be con- 


sidered as occupying the molecular orbital II,- 
But what is its contribution to the interaction po- 
tential? One might reasonably expect to find that 
since this single electron forms a one-electron bond, 
and since the strength of such a bond is usually 
about one-half that of the corresponding two- 
electron bond, a fairly good approximation ought 
to be 

re’ii)- j x a ] 2 j„ 

This question will now be more thoroughly in- 
vestigated in the light of the molecular orbital 
t heory. 

Using the ideas of Linnet t and allowing the 
extra electron to be located either in lb or Tb*, we 
have for the three-electron wave function 

O a (l)a(l) N a (l)ot(l) O a (l)/3(l) 

O a (2)«(2) N,(2)«(2) 0,(2)/3(2) 

0,(3)a(3) N a (3)a(3) 0 4 (3)/3(3) 

0 2 (l)a(l) N a (l)a(l) N a (l)0(l) 

±b O z (2)a{2) N a (2)a(2) N a (2)0(2) 

O a (3)a(3) N a (3)a(3) N z (3)j9(3) 

Wc expand these determinants and neglect the 
multiple exchange integrals giving rise to the pos- 
sibility of all three electrons being exchanged be- 
tween the two atoms. Such interaction terms arc 
usually quite small in comparison with other con- 
tributions, as can be sccti from overlap considera- 
tions. The interaction potential is then found 
to be 

U~— a 2 (0 2 N 2 :0JL ; 0/b)— 6 2 (0 2 X 2 :0 2 XyX 2 N 2 ) 

±2a6(0 < 0 I :N a N J :0 I N a ) 

where 

(0,X.:0,N,:0 i 0.)=J[. . . 0,(1)0.(2)0,(3)]* 

TI[. • - N I (l)N T J (2)0,(3)]rfr, etc. 

7/ is the Hamiltonian operator of the entire sys- 
tem, and the asterisk (*) signifies taking the com- 
plex conjugate. 

In general, if O and X are taken to represent 
symbolically any two atoms, there are two extreme 
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cases to consider as follows: 

(1) Suppose Unit X, for example, refers to an 
atom whose electron affinity is zero, so 
that the coefficient b vanishes. Normali- 
zation then requires that a be equal to 
unity and results in 

- (0 A AN ,:0 A) = - J Zi 

where the last equality follows from the 
definition of exchange. In other words, 
the quantity in parentheses refers to the 
situation in which two electrons are ex- 
changed back and forth between N and 
O, whereas the third electron remains sta- 
tionary on the 0 atom. Moreover, this 
interpretation is certainly consistent with 
the approximation of perfect pairing be- 
cause the spins are random, in which case 
we have 


f ~ 2 ( 2 <n)— Jzz 

(2) Suppose that N and 0 now refer to the 
same type of atom, for example oxygen, 
so that the coefficients, a and 6, arc each 
nearly equal to (2) -I/2 . The interaction 
potential then becomes 

V ~ — J„ ± (0,0, : N,N, AN,) 

where the last term signifies one electron 
on one atom, another electron on the 
other atom, and the third electron equally 
dividing its time between the two atoms. 
Thus, this last integral, together with its 
coefficient (unity in this case), is called a 
hybrid coulomb-ex change integral, and 
it clearly provides a measure of the im- 
portance of resonance configurations. 
Also, it is quite obvious from the form of 
the integrals involved that a good ap- 
proximation to V should be 

Y — — Jzzi^Jzz 

where X is a resonance parameter and, in 
general, depends on the interatomic sep- 
aration. Yanderslice, Mason, and Maisch 
(ref. 17 ) set. this quantity equal to a con- 
stant, 0.5, but it will be treated here as a 
variable. The plus sign, of course, refers 


to the “extra” electron being in the U> 
orbital and the negative to H z * occu- 
pation. 

Going hack to the original problem of the ground 
state of nitric oxide, it is certainly to be expected 
that the resonance parameter, X, will lie somewhere 
between its maximum value in case (2), in which 
the resonance is complete, and its value of zero 
in case (l). However, we still have to develop a 
procedure for its calculation. 

Since (OA:N,N,AN 2 ) must be greater than 
J zz because of overlap considerations, and since 
X must vanish whenever a or b is zero, a reasonable 
functional form seems to be 

X — • ya b 

where 7 in general depends on the interatomic 
separation and is greater than 2 in the case of 
complete resonance. 

To the extent that 7 can be regarded as constant, 
an assumption which will not he used in our future 
calculations, its value may be determined by a 
consideration of the problem described in case 
(2) above, in which a and b are both nearly equal 
to (2) -1/2 . Since this situation implies that the 
extra electron is just as likely to be found on O as 
on X, it can he imagined that the contribution of 
the three-electron bond to the interaction energy 
of the (n 2 ) 2 (n,*) state is 

that is, 7 — 3 (an excellent example is the resonance 
stabilization of IIe 2 + ). 

Proceeding one step further with this idea, it is 
now possible to derive an expression for the proba- 
bility (P) of the resonance ionic state 0 + — N~. 
This probability is equal to ¥ and, by using the 
normalization condition a 2 -\-b 2 =l, we can write 


where 


P~)-a 2 


_ A— _ x 

a ~'db~:iP ' /2 


Eliminating a and solving the ensiling quadratic 
equation yields 
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A somewhat better determination of the reso- 
nance parameter may be obtained by a considera- 
tion of the ground A^ 2 n state of nitric oxide in 
conjunction with the excited state, which will 
be described in the next section by use of the 
delta-function model. The electronic configura- 
tion of this latter state is 

so that the two energies involved arc 

T (TI) = Jxz^r^Jvv (29) 

and 

1 (lT^)^rx 2e Jyy (30) 

THE 6 y^ ; STATE OF NO 

A brief study of the potential function given in 
equation (30), especial!} 7 when cognizance is taken 
of the fact that J xx is usually many times greater 
than J vy in magnitude, clearly indicates that this 
particular state of nitric oxide is the most 
repulsive one obtainable from ground state atomic 
orbitals — -at least as long as the interatomic 
separation is not so small that the coefficients of 
linear combination of atomic orbitals are effec- 
tively different for Tfi and Tfi*. In addition, it 
should be noticed that <r type orbitals arc empha- 
sized, so that everything seems to point to this 
being the state most accurately described by a 
delta-function model with its associated liydrogen- 
like approximation of taking average electrons. 

As pointed out earlier, the situation in the case 
of hcteronuclear molecules is complicated some- 
what by the existence of two sets of delta-function 
strengths and shifting parameters. However, as 
will be justified later in connection with the treat- 
ment of exchange integrals, it is a good approxima- 
tion to take the geometric mean between corre- 
sponding homonuclcar molecules, 0 2 and N 2 in 
this problem, in which case the new potential 
function becomes 

F(‘S)= (n,K,)" 2 (27.206) g l9t ox p [- (f| * + c *»)] ev 

(31) 

where the subscripts 1 and 2 refer to oxygen and 
nitrogen, respectively; for example, 7q=6 and 
77 2 =5. 

The values of c h p i} and g t are, of course, given 


by expressions similar to those in equations (15) 
and, in particular, values for c 2 p 2 and g 2} as well as 
p, may be obt ained from the previous calculations 
oil nitrogen. Since the parameter p depends pri- 
marily on the number of electron pairs involved 
and seems to be fairly insensitive to a small 
change in this number, the assumption that 
p — 0.55 for 0 2 as well as X 2 is probably not too 
bad. Thus, the only additional information 
needed is the radius of the outermost electron shell 
of atomic oxygen, defined as the distance at which 
the electron charge density is a maximum and 
calculated to be 0.48 A by use of empirical screen- 
ing parameters (ref. 20) in Slater atomic orbitals, 
and the corresponding ionization potential (13.550 

°v). 

A straightforward calculation, using Newton's 
iteration procedure to solve equation (15c), 
yields the results presented in the third column of 
table V. Also shown in this table is the interac- 
tion energy corresponding to the A~ 2 IT state, as 
computed from the TIul 1 hi r t -Hi rschf elder function 
in equation (27). 

THE RESONANCE PARAMETER 

It will now be convenient to define a new param- 
eter k as the ratio of exchange integrals J xx to 
J V y. Equations (29) and (30) may then be 
written 

Vm=Jyy(k+\) (32) 

and 

VCJ:)=-J uu (k+2+k\) (33) 

Dividing the first by the second and solving for 
X yields 

, -\m*yz)+(k+ 2)TTn)i 
F( 6 XJ)+ifcF( i n) 

where everything on the right-hand side is known 
with the exception of k. Notice that X approaches 

TABLE V.— ENERGIES OF THE Pn AND f 'X! 

STATES OF NITRIC OXIDE 


T 

X 2 n 


2. 5 

-0. 1718 

0. 2093 

2. G 

1307 

. 2033 

2. 7 

- . 1009 

. 1533 

2. 8 

0788 

. 1150 

2. 9 

0021 

I 

. 0872 
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zero as k approaches infinity (sec table V). 

In order to determine this last parameter, use 
is made of the fact that the exchange energies 
should, in a major way, be dependent on the 
corresponding overlap integrals. In fact., a rea- 
sonable approximation would appear to be the 
following: 

J»= (CUsvO.NJ ~\sj[(O x O x :O x O z ) 

+ (N,N,:N,N.) +2(O x O x :N,NJ] (35) 

J„--- (0„X„:0„NJ »i S^WrOA) 

+ (N t N,:N,N,)+2(OA:N,N»)] (36) 

where 

(0*N*:0j.N r )= |t. • ■ Q,(1)0,(2)]* 

//[... N,(l)N x (2)Kr, etc. 

II is the complete Hamiltonian operator of the 
entire molecular system, and S u is the overlap 

integral jOj*N,(?r. 

Each of the integrals appearing in the right- 
hand members of equations (35) and (3G) repre- 
sents atomic plus coulombic energies. Since the 
coulombic energies may be shown to be essentially 
independent of whether p x or p y is used and to 
constitute only about 0.08 percent of the atomic 
energies for the distances under consideration, we 
have to a good approximation (since E 1 is small 


compared with E) 

(O x O x :O x O x ) « (0,0,10,0,) «2?( O) 

+E(N)+E l (0)-E l (K) (37) 

(X,N,:N,N*) - (N,N,:N,N,) «ff( O) 

+E(S)—E,(0) J rEi(S) (38) 

and 

(0*0* : N*N*) « (0,0, : N,N,) ~E(0)+E(K) (39) 

where E(0) and J?(N) are the total energies of 
atomic oxygen and nitrogen, respectively, and 
E x (0) and E X (K) refer to the corresponding one- 
electron energies. Thus, the expressions for J xx 
and J yy reduce to 



J„~S„*[E(0)+E(X)] 

(40) 

and 


J v »~S u i[E(0)+Em] 

(41) 

so that. 


/ e \2 

7 / \ 

* v ^ / 

\E>yy/ 

(42) 


The overlap integrals may be evaluated using 
elliptical coordinates in which the variables are 
(l/r)(r 0 +rx), v~ (1 //*) (r 0 — r N< ), and <p, and 
where r 0 and r N refer to the distances of an arbi- 
trary point from the nuclei O and N, respectively. 
The parameter r still represents the intemuelear 
separation and it is also the distance between the 
foci of the associated ellipse. 


Remembering now that / is in angstroms and using the ordinary atomic wave functions, we obtain the 
following expression for the parameter k: 


.1 

1X1 

p2s- 

[vV 2 — v 2 (m 4 +1)+M 2 ] exp (— j/i— yv) dip dv dp 

0 

^ 00 ^ i + 1 J 

(*2t 

1 

0 


2 — 1) — v 2 (p 4 — 1)+(m 4 — E)] exp (—.rp—yv) cos ?<p dtpdvdp 


(43) 


where 

x [(^ + ^)/4(0.5292)lr=3.9919r 

y [(>*! — 2 2 )/4 (0.5292)]/' =0.3 071/* 

z x effective nuclear charge seen by an L-shcll 
electron in atomic oxygen (—4.55) 


z 2 analogous quantity ( — 3.90) for atomic nitrogen 

In the case of fairly large intemuelear separa- 
tions, it is obvious from the definition of the 
elliptical coordinate system, and from the fact 
that the major contributions to the overlap 
integrals come from the region roughly half way 
between tlie nuclei, that a very reasonable 
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assumption at this point should be the setting of 
v equal to zero. For the purpose of obtaining an 
idea of the error thus introduced, we shall suppose 
that the maximum contribution is attained along 
the line (l/2)(a]+a 2 ) — a t «0.04 A to the oxygen 
side of the geometric center of the molecule. 
This corresponds to the position of touching of 
Bohr-type orbits, and the resulting shift will be 
assumed constant over all greater separations. 
Assuming, in addition, that the effective mag- 
nitude of ju remains of the order of unity, we have 
the following effects on the integrands of equation 
(43) : 

(1) The value of v for the above displacement 
of 0.04 A, and for r=2.7 A, is ^ —0.03, 
Since the smallest power of v appearing 
in the integrands is r 2 , it is easily seen 
that this particular effect of the assump- 
tion j/=0 is to increase only slightly the 
values of the overlap integrals. 

(2) The neglect of v , on the other hand, 
serves to decrease the values of the over- 
lap integrals through the omission of the 
factor exp(— yv) « exp ( + 0.025). 

(3) The exchange integrals are further de- 
creased when lull advantage is taken of 
the geometric mean procedure discussed 
in connection with equation (31). This 
implies that the expressions given in 
equations (40) and (41) should be multi- 
plied by the factor 

0 [/?(Q)/?(N)r 
E(0) + E(X) 

which is slightly less than unity. 

(4) Finally, the net effect of (1) and (2), 
only, is to decrease the exchange integrals 
by a very small amount for the range 
under consideration. However, as pointed 
out in a recent paper by the author (ref. 
7), and particularly in regard to the sup- 
posedly rigorous calculation of the A" 2 n 
state, the effective mi clear charges arc 
probably a bit larger than the ones men- 
tioned above. In consequence of this, 
the exchange integrals should be reduced 
slightly because of the more compact 
charge distributions and resulting de- 
crease in overlap. Thus, the elimination 
of v is perhaps better for our purposes 


than is a direct evaluation of the integrals 
in equation (43). 

In view of the above arguments, the expression 
for k may now be written 


k= 



fi 2 exp (— 


yt 2 ) exp (—x/z+Zm 


2 


or 


" (r 4 +2r 3 +2i’ 2 ) 1 2 

_(.r 3 +5j 2 + 12r + 12)_ 


(44) 


(45) 


As mentioned earlier, the character of the 
exponentials in equation (44), that is, the form of 
essentially a product of exp(— z x ?') and exp ( — z 2 r ) , 
justifies to a large extent the use of the geometric 
mean in the delta-function model in preference to 
the arithmetic average. 

The results for the resonance parameter, X, as 
calculated from equation (34) using equation 
(45), are presented in table VI and shown graph- 
ically in sketch (d). A comparison with the 
X=0.5 assumption of Vandcrslice, Mason, and 
Maisch (ref. 17) is also indicated. 


TABLE VI. THE RESONANCE PARAMETER FOR 
NO AND ASSOCIATED PROBABILITY 


r 

X 

P 

2. 5 

0. 5474 

0. 0345 

2. o 

. 5303 

. 0331 

2. 7 

. 5010 

. 0288 

2. 8 

. 4495 

. 0230 

2. 9 

. 3800 

. 0108 



Sketch (d). — Resonance parameter for NO. 
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Also shown in table VI is the probability of the 
resonance state 0 + — N~, as computed from equat ion 
(28). Its small magnitude further justifies our 
use of equations (32) and (33), without any ionic 
contributions, and explains the validity of the 
IIulburt-TIirschfelder curve fit in connection with 
the ground state. At the same time, since the 
IIulburt-TIirschfelder function also excludes any 
ionic effects, we can reasonably postulate that 
such good fortune will not prevail in the case of 
diatomic oxygen, where X is approximately equal 
to 1.5 and P is correspondingly much larger. 
These ideas have been confirmed in recent calcula- 
tions on 0 2 by the theoretical group at the 
University of Maryland (ref. 21). 

A question might reasonably be asked at this 
point as to why, if the ionic forces are negligible 
anyway, we even bother to introduce X into 
equations (32) and (33). The answer is, of 
course, in order to have some means of correcting 
or accounting for the fact, that the orbital occupied 
by the “extra” electron is warped one way or 
the other, depending upon whether it is an 
attractive or a repulsive bond, and thereby to 
include the perhaps much larger effect on the 
exchange forces. Thus, the use of a resonance 
parameter without an ionic force is entirely con- 
sistent in this problem. 

THE 4 n STATE OF NO 

As mentioned previously, Vanderslice, et id. 
(ref. 17), employ a very approximate method for 
obtaining the excited 4 II state of nitric oxide, with 
the result that a strange hump appears in the 
potential-energy curve even though no rotation is 
included. It was, therefore, thought interesting 
to pursue this question of the existence of an 
activation energy more thoroughly in the light 
of the present procedure. 

The electronic configuration of the l II state is 

(>*) 2 (Ib) (u*) (n,) 2 (n,*) 

with the corresponding energy given by 

U( 4 n) - J xx — (2 — X) J yy (46) 

or 


Dividing by equation (32), one readily obtains 

C( 4 n)=F( 2 n) (i-^) (48) 

The results are given in table VII, together 
with those of Vanderslice, et al., corresponding 
to points beyond the negative minimum, and it 
is seen that the bump has now disappeared. 

THE N 2 -0 2 INTERACTION 

The interaction between nitrogen and oxygen 
atoms, each of which is a member of its own 
respective diatomic molecule, lias already been 
discussed and summarized in equation (26), 
which may now be written 

F=-| J„(k+ 2) (49) 

Dividing by equation (32) and solving for V, 
we obtain 



the results of which are nicely curve fitted by the 
equation 

V=V 0 cxp(-ar) (51) 

The desired interaction between nitrogen arid 
oxygen molecules may now be expressed as a sum 
of the four atomic interactions, each of which is 
given in the form of equation (51), the orientation 
dependence being implicitly included through the 
four values of internudear separation. In many 
instances, however, and if the intermolecular 
separation is sufficiently large, it is convenient to 
average (ref. 17) the total interaction energy (4>) 

TABLE VII. THE V(r) INTERACTION FOR THE 
l n EXCITED STATE OF NITRIC OXIDE 


r 

Rtf. 17 

Eq. (48) 

2. 5 

0. 0488 

-0. 1657 

2. 6 

. 0453 

1265 

2. 7 

. 0376 

. 0979 

2. S 

. 0288 

. 0767 

2. 9 

. 0207 

~. 060G 


F( 4 n)- J w (ir— 2+X) 


(47) 
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over all orientations to obtain 
$(/*’)= 4 V 0 exp ( — a 7?) (a 3 /?(/ i</ 2 ) “ 1 

jh(a/?+2) sinli (“fi) sink (^) 
- 4 (? , !) s i" |, (^!)cos | i(^)] cv (52) 

where 

Jl distance in angstroms between the centers 
of mass of the molecules 
dijd -2 bond lengths of 0 2 and N 2 , respectively, that 
is, 1.207398 and 1.094 A (ref. 13) 

This equation yields the results shown in table 
VIII and plotted as curve III in sketch (e). Also 

TABLE VIII. -INTERACTION ENEEGIES FOR 
THE N 2 -0 2 SYSTEM 


R 

Eq. (52) 

Eq. (5.3) 

2 . 5 

0. 7377 

0. 7287 

2. G 

. 5054 

. 5G92 

2. 7 

. 4396 

. 4446 

2. 8 

. 3457 

. 3473 

2. 9 

. 274 1 

. 2713 



shown, for purposes of comparison, are the two 
computations of Vanderslice, Mason, and Maiseli 
(ref. 17), both with (curve II) and without (curve 
I) the addition of the second-order London dis- 
persion energy. The broken curve (V) is simply 
a graphically constructed connection between the 
present results and curve IV, which was derived 
(ref. 17) from measurements of the viscosities of 
diatomic nitrogen and oxygen at high tempera- 
tures. Finally, as illustrated in column 3 of table 
VIII, a fairly good approximation to equation (52) 
is given by the exponential formula 

<t>'(7?) =/l exp( — BI?) ev (53) 

where A and B have the values 350.3 and 2.470, 
respectively. 

As is evidenced by the somewhat greater con- 
sistency of curve III with the one derived from 
high-tempera ture-gas viscosity data and the rela- 
tively greater ease with which the gap can be 
covered by a reasonable interpolation, the present 
procedure of using the delta-function model, a 
variable resonance parameter, and the ratio of 
exchange integrals would seem to be preferable 
to that of previous calculations. A somewhat 
different analysis of the introduction of the pa- 
rameter k will be presented below in connection 
with 0 2 — 0 2 . 

The procedure of Vanderslice, Mason, and 
Maiseh, however, does present a reasonable first 
approximation and, in consequence of the fewer 
equations involved, may prove quite useful in 
dealing with more complicated species. In par- 
ticular, the choice of 0.5 for X is just about as good 
as can be obtained within the limitations imposed 
by the assumption that it is independent of inter- 
nuclear separation. The principal blame for the 
discrepancies between curves IT and ITT in sketch 
(e) probably lies in the fact that the range covered 
in the N 2 — 0 2 calculation falls precisely in the 
region of the hump in the 4 II state of nitric oxide, 
which state plays such a vital role in the Vaiulor- 
slice, et ah, method. 

It, is, on the other hand, difficult to assess the 
absolute accuracy of the present calculations be- 
cause of the lack of experimental scattering data 
to determine the potential curve at smaller sepa- 
rations. Its justification must depend to a large 
extent on the fact that the procedure involved is 
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basically the same as that of the previous treat- 
ment of N 2 — N 2j for which scattering data were 
available and with which excellent agreement was 
obtained. Indeed, the observable variations (e.g., 
the divergence of curve III from curve V resulting 
from the neglect of dispersion forces) may be ex- 
plained in a similar fashion to the explanation of 
the variations connected with the N 2 — X 2 problem. 
Moreover, a very interesting observation on curves 
I and II at 2.9 A indicates that their difference is 
very close to what is needed to make curve III 
coincide with V, that is, the simple London dis- 
persion energy may be added to the present results 
to obtain the complete potential curve beyond 
this point. 

OXYGEN INTERACTIONS 
METHOD 

The interaction energy between two oxygen 
molecules will now be determined in much the 


same manner as in the previous calculations on 
X 2 -X 2 and N 2 — 0 2 . Again, because of the 
compact charge distributions, 0 + — O - resonance 
(three-electron bonds) between atoms belonging 
to different 0 2 molecules will be ignored. There 
remains then the following set of nine possibilities: 



Oxygen atom 
(molecule A) 



Oxygen atom 
(molecule B) 


Number of electrons 
in— 

Number of electrons 
in — - 

p. 

Py 

p* 

p. 

Py 

p. 

P, 

Py 

p. 

P* 

1 V 

P t 

2 

1 

l 

1 

2 

1 

2 

1 

1 

i 

1 

2 

2 

1 

l 

1 

1 

2 

1 

2 

1 

2 

1 

1 

2 

1 

1 

1 

i 

2 

I 

1 

2 

1 

2 

1 

1 

1 

2 

2 

1 

1 

1 

i 

2 

; 

1 

2 

1 

1 

1 

I 

2 


A straightforward application of the approximation of perfect pairing to each of these configurations 
yields, since the spins are random, 



~ J z ^ + (V,*+ + (jxx+ g Jm-t + (jz x+ J»v+ 2 

+ Q Jxz+2Jm+ 2 + Q) J**) + (j** + \ + (jj 

+ Q J »+l J m + 2J„) 


(54) 


where the x direction is taken to coincide with the 
line of nuclear centers. 

On setting J vy —Ju as a consequence of axial 
symmetry, this equation may be reduced to 

— - {Jxz~\-~Jyv) (55) 

There remains now only the task of eliminating 
J xx and J V y from equation (55) in order to specify 
uniquely this interaction. As mentioned pre- 
viously, however, the ground A r3 S a “ state of 0 2 
cannot be used for this purpose because of the 
inability of a Hi dburt-IIirschf elder function to 
cope with the 0 + — 0~ resonance phenomenon 
arising from thrcc-clcctron bonds; that is, the use 
of this empirical function would yield a potential 
curve lying above the correct one. In fact, the 
only diatomic state, and there are IS (ref. 22) 


which can dissociate into normal QP) atoms, for 
which such an empirical function appears to be 
rigorously confirmed by a Rydbcrg-Klein-Recs 
calculation (spectroscopic data) is (ref. 21) 

3 A.:(^) 2 (n,) 2 (n,*) 2 (n,)(n/) 

where the two three-electron bonds of 

X 3 S i “I(0 2 (n,) 2 (n/)(n # ) 2 (n f *)] 

have been replaced by one consisting of two 
electrons and another containing four. 

The appropriate equations for this state are 

I 7 ( 3 AJ=0.9154[(1 — c -3 ) 2 

+ 0 .02 1 247£ 3 e - 2 ^( 1 + 1 .3282/5) - 1 ] (56) 

where 

(3=5.4037 (44q 4 4 ) (Hulburt-IIirschfclcIer) 
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and 

I C A u ) =J XX 3 Jyy 

(approximation of perfect pairing) (57) 

where the four-electron bond has been obtained 
by means of a straightforward extension of the 
procedure for the three-electron variety discussed 
previously in connection with nitric oxide. The 
result, of course, is that the electrons involved 
may again be described as belonging to atomic 
rather than molecular orbitals and with a con- 
tribution (random spins) of 4 (K — )tyJ yu =—2J vy 
to the interaction potential. The question of 
resonance thus does not enter this discussion. 

Finally, a delta-function computation per- 
formed on an excited 5 ll g state of 0 2 is used to 
complete the aforementioned 0 2 — 0 2 requirements 
in the long-range “tail” region. 

As a further check on the validity of the above 
representation of the 3 A U state, since there is some 
uncertainty in the numbering of the vibrational 
levels (ref. 21), a second calculation was performed 
using the 1 A g state as a basis. Despite the fact 
that the Rydberg-Klein-Rees curve for this state 
is not known over a large enough range to furnish 
a stringent test of the fit of an empirical potential, 
it is believed that the TIulburt-Hirschf elder func- 
tion should suffice because of the following two- 
and four-electron bond structures: 

'AA* x y(n v )\n*y(n 2 y 

The equations corresponding to (56) and (57) 
are 

F( 1 A,) = 4.230[(1 -e-*) 2 

+0.089501^-^(1 + 2.6976/3)-!] (58) 


where 


and 


0=3.4203 


>-1.2155\ 
k 1.2155 ) 


I 0+) — Jxx Jy\ 


(59) 


Altogether four determinations of T 7 (r) were 
made and they may be summarized as follows: 

1. 3 A U (ITulburt-Hirschfelder) and 5 U 0 (delta- 
function model) states of 0 2 

2. ! A„ (Hulburt-Hirschfelder) and 5 U 0 (delta- 
function model) states of 0 2 


3. 3 A U (Hulburt-Hirschfelder) and 1 A SJ (Hul- 
burt-Hirschfelder) states of 0 2 

4. 5 n o (delta-function model) state of 0 2 
and the ratio of exchange integrals tech- 
nique. An ulterior motive here is, of 
course, an investigation of the validity 
of the introduction of the parameter k in 
the N 2 — 0 2 problem. In addition, valu- 
able insight into the limitations of the 
approximation of perfect pairing should 
be gained. 

THE 5 U a STATE OF 0 2 

The electronic configuration most appropriate 
for treatment by a delta-function model is that 
state in which each valence electron is antibonding 
to its maximum extent, subject to the condition 
that the dissociation products are normal atoms. 
The state most closely fulfilling these requirements 
appears to be 

5 n 0 : (tr z ) (cr z *) 2 (n i/ ) (n„*) z (n 2 ) (n**) 

provided the internuclear separation is large 
enough that the magnitudes of t He coefficients of 
linear combination of atomic orbitals are not 
effectively different from 2“ 1/2 . 

In the case of complete resonance (homomidear 
molecules) it was found previously that the 
associated resonance parameter (X) should 
approach (for small separations) one of the values 
±3/2 for the “extra” electron in a three-electron 
bond. Actually, this statement is true over the 
entire range of r as long as one sta 3 ~s within the 
confines of simple molecular orbital theory; but 
second-order approximations, namely the intro- 
duction of configuration interaction, indicate that 
|X] should decrease from 3/2 at large distances. 
An investigation of these effects is presented in a 
later section. 

Tims, on remembering that a delta-function 
calculation does not in itself include ionic con- 
tributions, we obtain, by application of the 
approximation of perfect pairing, 

F( 5 n„)= (5«7„+7J w ) (60) 

The delta-function equations are basically the 
same as equations (15), except that the ionization 
potential (13.550 ev), the outermost atomic shell 
radius (0.48 A), and the effective number of 
hydrogen-like electron pairs (6) now refer to 
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oxygon rather than nitrogen. Again the param- 
eter p in equation (15b) is taken to be 0.55, which 
might possibly cause the interaction potential to 
be slightly too high; however, as explained above, 
this should be of only minor significance. The 
results of this computation, together with those 
from equations (56) and (58), are presented in 
table IX, 

THE V(r) INTERACTION 

Tf equations (55), (57), and (GO) are combined on 
the one hand, and (55), (59), and (60) on the other, 
the interaction between atoms belonging to two 
different 0 2 molecules may be expressed in the 
following two ways: 

(6i) 

and 

F(r)=| [2F( 5 n„)+F( , 4,)] ( 62 ) 

The results, as presented in the second and third 
columns of table X, are very consistent with each 
other, thereby indirectly implying that a Ilulburt- 
Hirsch folder function is valid for both t lie 3 A U 
and l A 0 states of 0 2 . 

The remaining two determinations of T T (r) are 
obtained from the combination of equations (55), 


TABLE IN. THE V(r) INTERACTIONS FOR 0-0 


r 

3 A u 

■A. 

5 IL 

2. 5 

- 0. 0387 

0. 0854 

0. 3278 

2. 0 

- . 0200 

0038 

. 2483 

2. 7 

0187 

0401 

. 1880 

2. 8 

0130 

0388 

. 1423 

2. 0 

0091 

0312 

. 1077 


TABLE X. -THE V(r) INTERACTION FOR 0 2 -0 2 


r 

Eq. (61) 

Eq. (02) 

Eq. (63) 

Eq. (04) 

2. 5 

0. 1278 

0. 1267 

0. 1382 

0. 1179 

2. G 

. 0071 

. 09G2 

. 1059 

. 0892 

2. 7 

. 0737 

. 0726 

. 0842 

. 0674 

2. 8 

. 0559 

. 0540 

. 0G8S 

. 0510 

2. 9 

. 0424 

. 0409 

. 0572 

. 038G 


(57), and (59), and from (55) and (60), together 
with the definition of k in equation (42), as 
follows : 

F(r)=^ (63) 

and 

FM=i“ F(X> (s£±f) «*> 


Since Zj = Z 2 in this case, the equation analogous to (43) becomes 




[fiV-^V+l)+M 2 l exp 

1 ) — ^ 2 (m 4 — 1 ) + (m 4 — M 2 )] oxp {—Xjj) cos 2 vd<pdvd)i 


V 


j 


(65) 


where z—Zrj 2(0.5292), and Z is the effective 
nuclear charge of an L-shell electron in atomic 
oxygen ( ^4.55), 

In contrast to the handling of the nitric oxide 
problem and because of the fact that y is 
identically equal to zero, thus making invalid one 
of the principal arguments in support of setting 
v equal to zero, equation (65) was evaluated 
directly to obtain 


/r 4 +2r*— 3 t 2 -15j-15Y 
\ / 3 +6r z +15r+15 / 


(06) 


Calculations based on equations (63) and (64) 
are shown as columns four and five of table X, 


where it is seen that they are in fairly good agree- 
ment with each other and with columns two and 
three. For example, the deviation between the 
results of equation (64) and the average of (61) 
and (62) is only about 7.5 percent, indicating that 
the ratio of exchange integrals technique is a 
reasonable approximation whenever sufficient spec- 
troscopic data are not available. A further 
discussion of this approach is presented in the 
section “Limitations of the Approximation of 
Perfect Pairing.” 

THE RESONANCE PHENOMENON 

Since the 5 n o state of (\. plays the role of a 
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common denominator in equations (61) and (62), 
and since the equations agree so well with each 
other and are at the same time at odds with the 
results of equation (63), the implications are quite 
strong that something must be wrong with either 
the delta-function model or the state chosen to be 
represented by it. 

An obvious first correction would seem to be a 
reduction in the parameter p of equation (15b); 
but it has been stated previously that, if anything, 
p should be greater than the X a value of 0.55 in 
order to approach the correct limit of unity as the 
atomic number increases. A more direct modifi- 
cation would be perhaps the addition of configura- 
tion interaction between states of the same sym- 
metry species to account for the breakdown of 
simple molecular orbital theory at large separa- 
tions. This effect is more important here than 
for N—N and N— O because of the more compact 
(larger effective nuclear charge) electronic charge 
distributions of O — O. Thus, the distances be- 
tween atoms need not be as great as before to 
produce the same molecular orbital deficiencies, 
a fact confirmed to some extent by the increase in 
percentage deviation from 8 to 27 between the 
results of equation (63) and the average of (61) 
and (62), as r goes from 2.5 to 2.9 A. 

Furthermore, the state, for example, 

described in equation (60) should interact to some 
degree with tin* configuration 

5 n u :(<rJVx*)(nj 2 (n/)(n,)(n/) 

the existence of which is to be expected on the 
basis of the incorrect simple molecular orbital 
prediction that both “extra” electrons might 
possibly be found on one atom, even at large 
separations. 

Although an exact treatment of this problem is 
beyond the scope of the present paper, it is never- 
theless apparent that the need for configuration 
interaction may be roughly satisfied by the intro- 
duction of a variable resonance parameter (j\|<^3/2 
and decreasing with increasing separation) into 
equation (60) ; that is, 

T (Tl 0 )= — (l-f-X)t/ xx (2-f-X)t/ yy (6/ ) 

where it is assumed that the X associated with the 
x direction is not far different from that associated 
with y. 

A simultaneous solution of this equation with 


(57) and (59) then yields 

F( 5 n„)-i [3C( 3 AJ-5TTA,)] 

x=_ rCAj-2F(‘A 0 ) (68) 

The results, as given in table XT and plotted in 
sketch (f), are of the predicted order of magnitude, 
and the general behavior is seen to resemble very 
closely that of the nitric oxide calculations of 
sketch (d). The variation in the latter, however, 
is smaller than that of equation (68) because, for 
example, of the existence of only one three-electron 
antibond in the f) ]C stale of XO. Consequently, 
conversion of this to a three-electron bond, as is 
done in the 0 2 ( 5 n„) state above, would further 
imply a change from approximately “gerade” to 
“ungcrade” symmetry. This results in smaller 
off-diagonal interaction matrix elements. 

Indeed, these agreements should prove quite 
influential in any future attempt to corroborate 
the initial assumptions. 

LIMITATIONS OF THE APPROXIMATION OF PERFECT PAIRING 

A determination of the parameter k from the 
3 A m and x \ states of 0 2 furnishes an excellent 

TABLE XI - THE 0 2 RESONANCE PARAMETER, X 


r 

Eq. (G8) 

2. 5 

1. 3049 

2. G 

I. 2828 

2. 7 

1. 1718 

2. 8 

1. 0044 

2. 9 

. 8147 



Sketch (f). — Resonance parameter for 0 2 . 
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example of the limitations of the approximation 
of perfect pairing. The simultaneous solution of 
equations (57) and (59), in conjunction with the 
definition of k, yields 


, F( 3 A w )-3F( J A ff ) 
F( 3 A a )-F( l A,) 


(69) 


the results of which are found in table XII, 
together with those from equation (G6), for 
purposes of comparison. 

The order of magnitude discrepancies, as well 
as the obviously incorrect decrease in k with 
increasing r, may be most easily explained in terms 
of the following three examples: 

(1) In the usual case of using the approxima- 
tion of perfect pairing to find the relation 
between three different states of a 
molecule, the equations involved arc of 
the general form 

I 1 = O.J rx + bJ yy-\-Fl 

T 2 — - CtJ XX T" (U \ yy “ {“ F 2 

T T 3 = eJ XX +jj yy + Fz 


where F i represents the corrections due 
to Weinbaum type (ref. 23) ionic terms, 
etc., in the wave function. It is then 
assumed that these corrections are small 
enough that they effectively cancel out 
in the equation of Vz y possibly, as a linear 
combination of V\ and V 2 ; or in other 
words, the corrections add in the same 
fashion as the exchange terms. The truth 
of these assertions is difficult to ascertain, 
although such seems to be the situation 
in the N 2 problem. On the other hand, 
this example may be too simple for any 
adequate judgment. 


(2) Another technique in which this approxi- 
mation might be employed is that of 
determining the relation between Fi, V 2 , 
and k , where the latter is obtained by 
means of an independent calculation. 
The appropriate equation in this case is 

V l ak+ b 4~ ( Fxjjyy ) Jlk+b 

V 2 (F2j jyy) ck-\-d 

provided k is large enough that the 
correction terms are insignificant in 
comparison. The agreement of the pre- 
vious N 2 — 0 2 curve with one obtained 
from high-tempera ture-gas viscosity data 
seems to justify this approach, at least 
insofar as states for which the ionic 
contribution is known to be of minor 
importance are concerned. 

(3) Finally, in the example of this section, 
k is found solely from the relation between 
Vi and V 2t but under the same general 
restrictions that apply in case (2). The 
ionic F functions associated with the 3 A U 
and l A ff states of O 2 , while apparently not 
strong enough to affect critically the 
Hulburt-Hirschf elder descriptions, never- 
theless are not so weak as to be negligible 
in the computation of k . In contrast to 
the ground state, where conversion of the 
two three-electron bonds to antibonds 
results in a quite different energy and 
ordinary empirical curve fits are not 
possible, t lie corresponding 3 A W and J A ff 
configuration interaction states have the 
same valence bond energy as their 
counterparts. These are, respectively, 

* a u : (^*) 2 (n,) 2 (n/)*(n I )(n 1 *) 

T J XX 3t7 yy 


TABLE NIL— THE RATIO OF EXCHANGE 
INTEGRALS, k 


r 

Eq (66) 

Eq. (G9) 

2. 5 

52. 56 

4. 655 

2. 6 

58. 71 

4. 455 

2. 7 

65. 22 

4. 227 

2. 8 

72. 10 

4. 009 

2. 9 

79. 33 

3. 820 


and 


(0 2 (n,) 2 (n/) 2 (ib*) 2 

M\Uy) 2 (u v *) 2 (u g ) 2 

(*.) 2 (n,)*( n,*m*) 2 

V— T — T 

» — wri yy 

Hence configuration interaction, together 
with its implied decrease in ionic forces, 
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is much more important for the latter 0 2 
excited species. 

In spite of these uncertainties, the 0 2 — 0 2 
interaction in the following section was computed 
on the basis of equation (63). A more valid 
representation, although not enough information 
is available to evaluate all the unknown param- 
eters, is 

V(r) (3TTAJ - 5F(‘A,) + 5F(>A f ) -SF(^u)] 

where the condition for the recovery of (63) is 
easily seen to be 


FQ A,)-!*’(*A„) 

The validity of our procedure must, therefore, rest 
on whether or not the three 1 A g configuration 
interaction states effectively lower the ionic 
contributions below those corresponding to Z A U . 
Such an occurrence does not seem too unreasonable 
in view of our present knowledge. 

In summary, the usefulness of the approxima- 
tion of perfect pairing lies in its ability to yield 
linear relations between potentials corresponding 
to several different states. Any other application, 
such as an attempt to analyze the component parts 
of a particular state or to calculate the energy 
directly, must be carried through only with the 
most extreme caution. 

THE 0 2 -0 2 INTERACTION 

Curve fitting the results of equation (63) by 
means of an exponential function, and then 
averaging the total 0 2 — 0 2 interaction energy over 
all orientations according to equation (20), where 
the London dispersion energy must be omitted 
and <7=1.207398 A for the bond length of ground 
state 0 2 , one obtains the potential values given in 
column 2 of table XIII. The results tabulated 
in column 3 were computed from the associated 
expression 


<£'(/?) = 146.6 exp (—2. 1097?) (70) 

and it is seen that the agreement here is^to about 
one significant figure. Finally, column 4 of this 
table represents the high-temperature-gas viscosity 
data (ref. 17) as calculated from the following 


TABLE XIII. — THE 0,-0* INTERACTION ENERGY 



modified Buckingham (exp-6) function (ref. 24) : 


$"(/?) 


€ 

1 — (6/«) 




(71) 


where 

e depth of the Van der Waals minimum (0.01137 
cv) 

I? m position of this minimum (3.726 A) 
a dimensionless parameter measuring the steep- 
ness of the repulsive energy (17.0) 

It is estimated from the temperature range of the 
viscosity experiments that this curve is not valid 
for distances smaller than about 3.2 A. 

In sketch (g) are plotted the results of the pres- 
ent paper (I), together with those from the vis- 
cosity measurements (II), and the ease with which 
the gap between the two regions can be covered by 
a reasonable interpolation (III) is quite apparent. 
The divergence between curves I arid ITI in the 
neighborhood of 2.9 A may again be attributed to 



Sketch (g). — Comparisons of 0 2 — 0 2 interaction energies. 
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the neglect of dispersion contributions, effects 
which are more important here than in the pre- 
vious examples because of the smaller cutoff dis- 
tance arising from the interpenetration of more 
closely packed charge distributions. 

Also shown in sketch (g) is a curve (TV) derived 
from vibrational relaxation time data (refs. 25—29), 
this being the only experimental information avail- 
able in the region covered by the present theory. 
An analysis of vibrational deactivation theory 
shows that the discrepancies between curves I and 
IV are to be expected on the basis that the latter 
is a one-dimensional treatment involving only the 
end-to-end molecular configuration. This geome- 
try is chosen, of course, because of the maximum 
energy transfer on collision from vibration to trans- 
lation; however, it also represents the largest inter- 
action energy of any molecular orientation and 
should not strictly be compared with t he random 
situation covered here. An additional basic defect 
of the vibrational theory, as pointed out by Cot- 
trell and Ream (ref. 30), is that the slope of the 
true interaction curve is not nearly as large as that 
of the theoretical one, an observation which is 
clearly borne out in this example. 


The recent calculations of Vanderslice, Mason* 
and Maisch (ref. 21) have been omitted from tbe 
above comparisons because of essentially the same 
reasons discussed in the preceding paragraph. In 
effect then, any absolute confirmation of the be- 
havior of curve I must await the performance of 
a suitable scattering experiment. 

A final rather interesting observation is that the 
N 2 — N 2 repulsive curve is generally lower than the 
one corresponding to 0 2 — 0 2 . Thus, the effect of 
the existence of one more electron pair in any 
O—O type of interaction seems to overshadow the 
fact that the electronic charge distribution of nitro- 
gen is less compact and thereby results in greater 
individual repulsive overlap. The N 2 — 0 2 curve, 
on the other hand, tends to lie slightly lower than 
either of the above in the short-range region 
(<2.6 A) because tbe overlap between similar 
transverse orbitals A and B on atoms of different 
species is generally smaller than that between 
orbitals A- A or B-B at the same separation. 

Ames Research Center, 

National Aeronautics and Space Administration, 
Moffett Field, Calif., Sept. 25, 1059. 
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